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ABSTRACT 

We study the tidal excitation of gravity modes (g-modes) in compact white dwarf bi- 
nary systems with periods ranging from minutes to hours. As the orbit of the system 
decays via gravitational radiation, the orbital frequency increases and sweeps through 
a series of resonances with the g-modes of the white dwarf. At each resonance, the 
tidal force excites the g-mode to a relatively large amplitude, transferring the orbital 
energy to the stellar oscillation. We calculate the eigenfrequencies of g-modes and their 
coupling coefficients with the tidal field for realistic non-rotating white dwarf models. 
Using these mode properties, we numerically compute the excited mode amplitude in 
the linear approximation as the orbit passes though the resonance, including the back- 
reaction of the mode on the orbit. We also derive analytical estimates for the mode 
amplitude and the duration of the resonance, which accurately reproduce our numer- 
ical results for most binary parameters. We find that the g-modes can be excited to a 
dimcnsionlcss (mass- weighted) amplitude up to 0.1, with the mode energy approach- 
ing 10~ 3 of the gravitational binding energy of the star. Therefore the low-frequency 
(<; 10~ 2 Hz) gravitational waveforms produced by the binaries, detectable by LISA, 
are strongly affected by the tidal resonances. Our results also suggest that thousands 
of years prior to the binary merger, the white dwarf may be heated up significantly 
by tidal interactions. However, more study is needed since the physical amplitudes of 
the excited oscillation modes become highly nonlinear in the outer layer of the star, 
which can reduce the mode amplitude attained by tidal excitation. 
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1 INTRODUCTION 

It is well known that non-radial gravity modes (g-modes) are 
responsible for the luminosity variations observed in some 
isolated white dwarfs (called ZZ Ceti stars) in the instabil- 
ity strip. These g-modes are thought to be excited by a con- 
vective driving mechanism operating in the shallow surface 
convection zone of the star (see Brickhill 1983; Goldreich & 
Wu 1999; Wu & Goldreich 1999). 

In this paper we study the tidal excitation of g-modes in 
compact binary systems containing a white dwarf (WD) and 
another compact object (white dwarf, neutron star or black 
hole). The Galaxy is populated with ~ 10 8 WD- WD bina- 
ries and several 10 6 of double WD-NS binaries [Nelemans et 
al. (2001); see also Nelemans (2009) and references therein]. 
A sizeable fraction of these binaries are compact enough so 
that the binary orbit will decay within a Hubble time to 
initiate mass transfer or a binary merger. Depending on the 
details of the mass transfer process (including the response 
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of the WD to mass transfer), these ultra-compact binaries 
(with orbital period less than an hour) may survive mass 
transfer for a long time or merge shortly after mass trans- 
fer begins. A number of ultra-compact interacting WD- WD 
binary systems have already been observed [including RX 
J0806. 3+1527 (period 5.4 min) and V407 Vul (period 9.5 
min); see Strohmayer 2005 and Ramsay et al. 2005]. Re- 
cent surveys (e.g., SDSS) have also begun to uncover non- 
interacting compact WD binaries (e.g., Badenes et al. 2009; 
Mullally et al. 2009; Kilic et al. 2009; Marsh et al. 2010; 
Kulkarni & van Kerkwijk 2010; Steinfadt et al. 2010). De- 
pending on the total mass, the systems may evolve into Type 
la supernovae (for high mass), or become AM CVn binaries 
or R CrB stars (for low mass). Many of these WD binaries 
are detectable in gravitational waves by the Laser Interfer- 
ometer Space Antenna (LISA) (Nelemans 2009). 

In this paper we consider resonant tidal interaction in 
WD binaries that are not undergoing mass transfer. This 
means that the binary separation D is greater than D m i n , the 
orbital radius at which dynamical merger or mass transfer 
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occurs, i.e. 
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where M is the WD mass, Mt — M + M' is the total mass 
and R is the WD radius. This corresponds to orbital periods 
of 
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Since WD g-mode periods are of order one minute or longer, 
they can be excited by the binary companion prior to 
mass transfer. In particular, as the binary orbit decays 
due to gravitational radiation, the orbital frequency sweeps 
through a series of g-mode frequencies, transferring orbital 
energy to the modes. Although the overlap integral of the 
g-mode eigenfunctions with the tidal potential is generally 
quite small, a binary system that spends a long time at res- 
onance can still excite g-modes to large amplitudes. 

Previous studies of tidal interaction in WD binaries 
have focused on quasi-static tides (e.g. Iben, Tutukov & 
Fedorova 1998; Willems, Deloye & Kalogera 2009), which 
essentially correpond to non-resonant f-modes of the star. 
Such static tides become important only as the binary ap- 
proaches the tidal limit [equation Q]. Racine, Phinney & 
Arras (2007) recently studied non-dissipative tidal synchro- 
nization due to Rossby waves in accreting ultra-compact 
WD binaries. Rathore, Blandford & Broderick (2005) stud- 
ied resonant mode excitations of WD modes in eccentric 
binaries. They focused on f-modes, for which the resonance 
occurs when harmonics of the orbital frequency matches the 
mode frequency. As mentioned above, for circular orbits, 
such resonance with the f-mode does not occur prior to mass 
transfer or tidal disruption. Their published analysis also did 
not include back reaction of the excited mode on the binary 
orbit. 

The problem of resonant mode excitations in compact 
binaries has been studied before in the context of coalescing 
neutron star binaries: Reisenegger & Goldreich (1994), Lai 
(1994) and Shibata (1994) focused on the excitations of g- 
modes of non-rotating neutron stars; Ho & Lai (1999) and 
Lai & Wu (2006) studied the effects of NS rotation - includ- 
ing r-modes and other inertial modes; Flanagan & Racine 
(2006) examined gravitomagnetic excitation of r-modes. In 
the case of neutron star binaries, the orbital decay rate (for 
orbital frequencies larger than 5 Hz) is large and the mode 
amplitude is rather small, so the back reaction of the excited 
mode on the orbit can be safely neglected (see section 5 of 
the present paper). By contrast, in the case of WD binaries, 
the orbital decay is much slower and the excited mode can 
reach a much larger amplitude. It thus becomes essential to 
take the back reaction into account. 

In this paper, we consider WD binaries in circular or- 
bits, consistent with the observed population of compact 
WD binaries (e.g. Kulkarni & van Kerkwijk 2010). Such cir- 
cular orbits are a direct consequence of the circularization by 
gravitational radiation and/or the common envelope phase 
leading to their formation. A key assumption of this pa- 
per is that we assume the WD is not synchronized with 
the binary orbit. While it is true that the tidal circular- 
ization time scale is much longer than the synchronization 
time, the observed circular orbit of the WD binaries does 



not imply synchronization. While there have been numer- 
ous studies of tidal dissipation in normal stars and giant 
planets (e.g., Zahn 1970,1989; Goldreich & Nicholson 1977; 
Goodman & Oh 1997; Goodman & Dickson 1998; Ogilvie & 
Lin 2004,2007; Wu 2005; Goodman & Lackner 2009), there 
has been no satisfactory study on tidal dissipation in WDs. 
Even for normal stars, the problem is not solved (especially 
for solar-type stars; see Goodman & Dickson 1999; see Zahn 
2008 for review). In fact it is likely that the excitations of 
g-modes and other low-frequency modes play a role in the 
synchronization process. The orbital decay time scale near g- 
mode resonances is relatively short (of order 10 4 years for or- 
bital periods of interest, i.e., minutes), so it is not clear that 
tidal synchronization can compete with the orbital decay 
rate. Given this uncertainty, we will consider non-rotating 
WDs (or slowly-rotating WDs, so that the g-mode properties 
are not significantly modified by rotation) as a first step, and 
leaving the study of the rotational effects to a future paper. 

This remainder of the paper is organized as follows. In 
section 2 we present the equations governing the evolution of 
the orbit and the g-modes. Section 3 examines the properties 
of WD g-modes and their coupling with the tidal gravita- 
tional field of the companion. In section 4, we numerically 
study the evolution of the g-modes through resonances, and 
in section 5 we present analytical estimates of the resonant 
g-mode excitation. We study the effect of mode damping on 
the tidal excitation in section 6 and discuss the uncertainties 
and implications of our results in section 7. 



2 COMBINED EVOLUTION EQUATIONS FOR 
OSCILLATION MODES AND BINARY 
ORBIT 

We consider a WD of mass M and radius R in orbit with 
a companion of mass M' (another WD, or NS or BH). The 
WD is non-spinning. The gravitational potential produced 
by M' can be written as 
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where r = (r, 9, <j>) is the position vector (in spherical coordi- 
nates) of a fluid element in star M, D(t) = (D(t),n/2, $(i)) 
is the position vector of M' relative to M (D is the binary 
separation, $ is the orbital phase or the true anomaly) and 
the coefficient Wi m is given by 
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(Here the symbol (— ) p is zero if p is not an integer.) The 
dominant / = 2 tidal potential has W2±2 = (3n /lO) 1 ^ 2 , 
W20 = (tt/5) 1 ' 2 , W%±i — 0, and so only the m = ±2 modes 
can be resonantly excited. 

The linear perturbation of the tidal potential on M is 
specified by the Lagrangian displacement £(r,t), which sat- 
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isfies the equation of motion 
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where C is an operator that specifies the internal restoring 
forces of the star. The normal oscillation modes of the star 
satisfy £• £ Q = oj 2 £ a , where a — {n, I, m} is the usual mode 
index and io a is the mode frequency. We write £(r, t) as the 
sum of the normal modes: 
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The (complex) mode amplitude a a (t) satisfies the equation 
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where Q a is the tidal coupling coefficient (also used by Press 
& Teukolski 1977), defined by 
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Here Sp a = — V ■ (p£ a ) is the Eulerian density perturbation. 
In deriving Q we have used the normalization 
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Resonant excitation of a mode a occurs when uj a = mtt, 
where fi is the orbital frequency. 

In the absence of tidal interaction/resonance, the WD 
binary orbit decays due to gravitational radiation, with time 
scale given by (Peters 1964) 
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where Mt = M + M' is the total binary mass. When a 
strong tidal resonance occurs, the orbital decay rate can be 
modified, and we need to follow the evolution of the orbit 
and the mode amplitudes simultaneously. The gravitational 
interaction energy between M' and the modes in star M is 
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where we have restricted to the I = 2 terms and set G = 1. 
The orbital evolution equations, including the effects of the 
modes, are then given by 
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The last terms on the right-hand side of equations ( 12 1 and 



( 13 1 are the leading-order gravitational radiation reaction 
forces, with (see Lai & Wiseman 1996 and references therein) 
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where p = MM'/M t and v 2 = D 2 + (L»$) 2 . In equations 



( 12 H 15 1 we have set G = c = 1. We have dropped the other 
post-Newtonian terms since they have negligible effects on 
tidal excitations. The mode amplitude equation is given by 
equation |7f, or, 

M'W lm Q a 
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3 WHITE DWARF G-MODES AND TIDAL 
COUPLING COEFFICIENTS 

The non-radial adiabatic modes of a WD can be found by 
solving the standard stellar oscillation equations, as given 
in, e.g., Unno et al. (1989). The g-mode propagation zone 
in the star is determined by uj 2 , < iV 2 amd w 2 < L 2 , where 
Li = + Tja s /r is the Lamb frequency (a s is the sound 
speed), and iV is the Brunt- Vaisala, frequency, as given by 
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where g the gravitational acceleration, and the subscript "s" 
means that the adiabatic derivative is taken. Alternatively, 
iV 2 can be obtained from (Brassard et al. 1991) 
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The Ledoux term B accounts for the buoyancy arising from 
composition gradient: 

X y dlnY 
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\t dlnP ' 



where 



Xy 



I ainF J 



(21) 



(22) 



P,T 



and Y is the mass fraction of helium. This equation is valid 
for a compositional transition zone containing helium and 
one other element, as is the case for typical compositionally 
stratified DA WD models. 

Figure [T] shows the profiles of the Brunt Vaisala and 
Lamb frequencies for one of the WD models adopted in 
this paper. These models were provided by G. Fontaine (see 
Brassard 1991). Since the pressure in the WD core is al- 
most completely determined by electron degeneracy pres- 
sure, iV 2 oc xt is very small except in the non-degenerate 
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Table 2. Same as table 1, for a WD model of identical mass and 
composition but with T fj = 5080A". 



0.4 0.6 

Radius (r/R) 

Figure 1. The square of the Brunt Vaisala (solid line) and Lamb 
(dotted line) frequencies and the density (thick solid line) as a 
function of normalized radius in a DA WD model, with M = 
O.6M , R = 8.97 X 10 3 km, T cff = 8720 K. The spikes in the 
Brunt Vaisala frequency are caused by the composition changes 
from carbon to helium, and from helium to hydrogen, respectively. 



Table 1. The eigenfrequency ui a , tidal overlap parameter Q a , 
and numerical f-mode overlap co for the first six 1 = 2 g-modes of 
a white dwarf model. The white dwarf model has T E g = 8720A", 
M = 0.6M Q , and R = 8.97 X 10 3 km. Note that Q a and Q a 
are in dimensionless units such that G = M = R = 1, and 
(GM/i? 3 ) 1 / 2 /(2 7 r) = 0.053 Hz. 
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outer layers. As a result, g-modes are confined to the outer 
layers of the star below the convection zone. Lower-order 
modes have higher eigenfrequencies, so they are confined to 
regions where N 2 is especially large, i.e., just below the con- 
vection zone. Higher order modes have lower eigenfrequen- 
cies and can thus penetrate into deeper layers of the star 
where the value of A 2 is smaller. Cooler WDs have deeper 
convection zones that cause the modes to be confined to 
deeper layers where N 2 is smaller. Consequently, the eigen- 
frequencies and associated values of Q a tend to be smaller 
in cooler WDs due to the decreased value of A" 2 in the region 
of mode propagation. 

The other feature of WDs that strongly effects their g- 
modes is their compositionally stratified layers. The sharp 
composition gradients that occur at the carbon-helium tran- 
sition and the helium-hydrogen transition create large val- 



ues of the Ledoux term B [equation (21 1], resulting in sharp 
peaks in N 2 as seen in Figure [l] These peaks have a large 
effect on the WD g-modes, leading to phenomena such as 
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mode-trapping (e.g., Brassard 1991) and irregular period 
spectra. Thus, the eigenfrequencies and eigenfunctions of 
WD g-modes are very sensitive to WD models. 

Tables 1-2 give the / = 2 f-mode and g-mode frequen- 
cies and their tidal coupling coefficients for two WD models. 
While the full oscillation equations need to be solved to ac- 
curately determine the f-modes, the Cowling approximation 
(in which the perturbation in the gravitational potential is 
neglected) gives accurate results for g-modes. Since high- 
order g-modes have rather small \Q a \, the mode eigenfunc- 
tion must be solved accurately to obtain reliable Q a - To 
ensure that this is achieved in our numerical integration, we 
use the orthogonality of the eigenfunctions to check the ac- 
curacy of the value of Q a (see Reisengger 1994 for a study 
on the general property of Q a ). Since the numerical deter- 
mination of an eigenfunction is not perfect, it will contain 
traces of the other eigenfunctions, i.e., 
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with c a ~ 1 and \cp\ <C 1 for /? / a. This means that the 
numerical tidal overlap integral is 
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Since \Qo\ (for the f-mode) is of order unity, while \Q a \ <C 1 
for g-modes, to ensure (Qar)num accurately represents the 
actual Q a , we require 
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The results shown in tables 1-2 reveal that |co| is always 
more than an order of magnitude less than Q a , so the above 
condition is satisfied for the modes computed in this paper. 

We note from tables 1-2 that while in general higher- 
order g-modes tend to have smaller \Q a \, the dependence of 
\Q a \ on the mode index n is not exactly monotonic. This 
is the result of the mode trapping phenomenon associated 
with composition discontinuities in the WD. To see this, we 
note that a mode with amplitude £ a has energy given by 
E a = to a f d 3 x p |£ Q | 2 , thus we can define the mode energy 
weight function 

dE a 



<pr £ + J(i + l)G. ]H V 



(26) 



dlnP 

where H v = dr/dlnP = P/(pg) is the pressure scale height, 
and we have used 
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Figure 2. The mode energy weight functions for the n = 1 
(thickest line), n = 2 (thick line), and n = 5 (top line) modes 
(all for I = 2) for a WD with T cff = 8720K, M = 0.6M e ,R = 
8.97 X 10 3 km, displayed as a function of logP so that the structure 
of the outer layers of the WD is more evident. The y-axis for 
a given mode is intended only to show the relative value of the 
weight function. The squares of the Brant Vaisala (thin solid line) 
and Lamb (dotted line) frequencies are displayed to demonstrate 
how their values constrain the region of mode propagation. 



(e r is the unit vector in the r-direction). Figures [2] and [3] dis- 
play the weight functions for several g-modes of WD models. 
We can see that the weight functions for all the low-order 
modes are largest in the region below the convective zone 
near the spikes in TV 2 produced by the composition gradi- 
ents. For the modes shown in Figure [2] the smooth fall-off 
of the weight function just below the convective zone indi- 
cates that these modes are confined by the falling value of 
the Lamb frequency in this region. The weight functions of 
higher-order modes and the modes in WDs with deeper con- 
vective zones may drop sharply at the convective boundary, 
indicating that these modes are trapped by the convective 
zone rather than the decreasing Lamb frequency. 

The weight functions also reveal the phenomenon 
known as mode trapping caused by the composition gra- 
dients. Mode trapping is especially evident for the n = 2 
mode, as it is confined to the helium layer between the 
spikes in iV 2 . It is clear that the mode is reflected by the 
carbon-helium boundary at larger depths and by the helium- 
hydrogen boundary at shallower depths. See Brassard (1991) 
for a more detailed description of the effects of mode trap- 
ping. 

The weight function is essentially the energy of a mode 
as a function of radius, so it tells us where orbital energy is 
deposited when a mode is excited. Since the weight function 
is largest in the hydrogen and helium layers just below the 
convection zone, most of the mode energy exists in this re- 
gion of the WD. Thus, if the mode is damped, most of the 
mode energy will be damped out in this region. 
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Figure 3. Same as Figure [2] except for a WD model with 
Tcff = 5080K. Note that the convection zone extends deeper in 
this model, pushing the modes to larger depths. 



4 NUMERICAL RESULTS FOR MODE-ORBIT 
EVOLUTION THROUGH RESONANCE 

Having obtained the mode frequency and the tidal coupling 
coefficient, we can determine the combined evolution of the 
resonant mode and the binary orbit using equations Q, 
(12 1, and (131. These are integrated from well before reso- 



nance until well after the resonance is complete. The initial 
mode amplitude b a and its derivative b a (prior to a reso- 
nance) are obtained by dropping the b a , b a and fl terms in 
equation (16 1, giving 
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with Cl ~ — (3D/2D)Q. These expressions are valid for for 
(uJa/mQ) 2 - 1 > tl/(mQ, 2 ) ~ 3/(2mf2t_o) (see section 5). 

The evolution equations Q, ( |12[ ), and ( |13[ ) form a very 
stiff set of differential equations. The reason for this is that 
the problem involves two vastly different time scales: the or- 
bital decay time scale which is on the order of thousands 
of years, and the orbital time scale (or the resonant mode 
oscillation period) which is on the order of minutes. Con- 
sequently, a typical Runge-Kutta scheme would require the 
integration of millions of orbits, demanding a high degree of 
accuracy for each orbit. To avoid this problem, we employ 
the Rossenberg stiff equation technique (Press et al. 2007). 
The integrator requires a Jacobian matrix of second deriva- 
tives, meaning that we need to supply the 8x8 matrix of sec- 
ond derivatives corresponding to our 8 first-order differential 
equations. The evolution equations are sufficiently simple 
that this matrix can be found analytically. 

Figure [4] depicts an example of the mode amplitude and 
orbit evolution near the resonance for the n = 3 g-mode. Be- 
fore the resonance, the mode must oscillate with the same 
frequency as the binary companion, so the amplitude \b a \ 
is smooth. After the resonance, the mode oscillates at its 
eigenfrequency (which is now different from the forcing fre- 
quency). The amplitude of the mode continues to fluctuate 
after resonance because it is still being forced by the bi- 
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Figure 4. Mode amplitude |b a | and orbital distance D as a 
function of time during resonance. The amplitude and distance 
oscillate after resonance, causing the curves to appear as filled 
shapes due to the short period of the oscillations with respect 
to the total integration time. These oscillations occur due to the 
continued interaction between the excited mode and the orbit. 
Note the sharp drop in orbital distance at resonance, which is 
caused by the transfer of orbital energy into the mode. The mode 
parameters are given in table I (n = 3). The companion mass is 
M' = M = O.6M . 



Figure 6. The average post-resonance amplitude \b a \ for the 
n = 4 mode given in Table I as a function of the mass of the 
binary companion. The post-resonance amplitude increases with 
the mass of the binary companion as predicted by the analytical 
estimate except for very high companion masses. 

with an abnormally high value of Q a (trapped modes) may 
reach large amplitudes as well. 



< 



Figure 5. The average post-resonance amplitudes for the first 
five modes given in Table I. The open triangles mark the results 
obtained from numerical integration, and the asterisks indicate 
the results predicted by the analytical estimate described in sec- 
tion [5] The analytical estimates are usually accurate within a 
factor of 10%. The n = 2 and n = 5 modes are trapped modes. 



nary companion, although the amplitude of these fluctua- 
tions diminishes over time as the orbit moves further from 
resonance. 

Figure[5]displays the average post-resonance amplitudes 
for the first five modes given in Table I. Note that no mode 
exceeds a maximum amplitude of \b a \ = 0.1, and we expect 
that our linear approximation is a reasonable first approach 
to the problem before non-linear effects can be included. In 
general, lower-order modes reach larger amplitudes (due to 
their larger coupling coefficients), but higher-order modes 



5 ANALYTIC ESTIMATE OF THE 
RESONANT MODE AMPLITUDE 

Here we provide an analytic estimate of the mode amplitude 
attained during a resonance as well as the temporal dura- 
tion of the resonance (i.e., the characteristic time during 
which the resonant mode receives most of its energy from 
the orbit). 

For a WD oscillation mode with frequency ui a , the res- 
onant orbital radius is 

Prior to the resonance, as the orbital radius D decreases, 
the mode amplitude grows gradually according to equation 
(28 1. At the same time, the orbit also loses its energy to 
gravitational waves (GWs) at the rate 



Egw = 



32(MM') 2 M t MM'/2D 



5c 5 D 5 



to 



(31) 



where to is given by equation < |10j . We can define the be- 
ginning of resonance as the point where the orbital en- 
ergy is transferred to the mode faster than it is radiated 
away by GWs. That is, the resonance begins at the radius 
D — D a +(> D a ), as determined by 



E a = LBg'V 



(32) 



where E a is the energy contained in the mode. Near the reso- 
nance, the mode oscillates at the frequency close to ui a , so we 
can write the mode energy as E a ~ 2<J^b 2 a (including both 
the m — 2 and m = —2 terms), assuming uj a ^> \bg /b a \- 
Thu s we have E a ~ Auj\b a b a . Using equations (281 and 
( 29 1 for b a and b a , we find 



oj^ — (mil) 2 = 



2AujiM'{W lm Q a f 



MD 5 



1/3 



at D = D 



(33) 
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D a 



1 + 0.436 



ML 

~M 



1/3 



(Mi 



-5/9 



-4/9^2/3 



(34) 



where we have used I = m = 2 and ui a , Q a are in dimen- 
sionless units where G — M = R — 1. The mode energy at 
D = D a + is 



E a (D a 



=0.0701 



-8/9 



„ -10/9 
X UJ, 



(W lm Q a ) 



2/3 



M 

IT 



(35) 



Since for D < D a +, the orbital energy will be deposited 
into the mode much faster than it is being radiated away, 
we approximate that all the orbital energy between D a + 
and D a is transferred to the mode. Thus the mode energy 
increases by the amount 

(MM' MM'' 



AE a = 2 x 



V 2D a 2D 

Ot-\ 



(36) 



where we have multiplied by a factor of two to account for 
the fact that energy is also deposited at a nearly equal rate 
before resonance as it is after resonance. Using equation 
J34l, we find 



AE a =0.2804 



M'\ 4/3 ( M t 



M 



M 



-8/9 



R 



(37) 



This is exactly four times of equation ( 35 1. Thus the maxi- 
mum mode energy after resonance is E a ,max — E a {D a -\-) ~t~ 
AE a , or 



E, 



a,max 



5.75 x 10"" 



M^_\ i/3 fM 
M ) \M 



1 0.2 J \W- 3 J 



2/3 



M 2 

w 



(38) 



The corresponding maximum mode amplitude is 

-4/9 

&«,max =±8.48 • 10 ' " ' 



2 ,'Af\ 2/3 /Mt 
M ) \M 



0.2/ U0- £ 



1/3 



(39) 



Figures [5] and [6] compare our numerical results with the 
analytical expressions (j38h- ( 39 1 . We find good agreement for 



all the WD resonant modes considered. Figure [5] verifies the 
dependence of & a , m ax on the mode frequency and the value 
of Q a , while figure [fi] verifies the dependence of b a ,m&x on 
the mass of the binary companion (except for the highest 
mass cases discussed below). Therefore equations ( |38[ ) and 
( 39 1 provide fairly accurate estimates of the mode amplitude 
and energy without performing numerical integrations. 

For the very high companion masses (M' <; 10 3 M) 
shown in Figure [6] our analytical formula significantly over- 
estimates the post-resonance amplitude. The reason for this 
is that the gravitational decay time scale is shorter if the 
companion is more massive. If the companion is massive 
enough, the orbit will decay through resonance due to grav- 
itational radiation before the orbital energy of equation ( |37[ ) 
can be deposited in the mode (see below). Consequently, the 



amplitude to which a mode is excited decreases if the mass 
of the companion becomes very high. Therefore our analyt- 
ical formula overestimates the post-resonance amplitude for 
extremely massive companions. For any reasonable WD or 
NS masses, our analytical estimate is accurate, but for a 
super-massive black hole the estimate may become inaccu- 
rate. 

It is interesting that the above analytical results for the 
resonant mode energy is independent of the gravitational 
wave damping time scale to, in contrast to the NS/NS or 
NS/BH binary cases. In fact, for the above derivation to be 
valid, the following four conditions must be satisfied at D a +- 



(i) oj q > |& Q /6 Q | , 

(ii) uj a — (mfi) 2 ^> mSl, 

(iii) ui a — (mO,) 2 ^> 2mQ. \b a /b a 

(iv) uia — (mQ.) 2 » 6a/&J . 
Conditions (i) and (ii) both lead to 

u a tD 3> 3 
condition (iii) gives 

LO a to 3> 6 

and condition (iv) yields 

0J a t D > VT8 



{mQ,y 



{mQ,) 2 



(mfi) 2 



3/2 



(40) 
(41) 
(42) 

(43) 
(44) 

(45) 

(46) 



Clearly, condition (iii) is most constraining. 



In equations (44 1- ( 46 1, the right-hand sides should be eval- 
uated at D a 
With 

Ji - (mf2) 2 



0.0064 



= 3 


(Da 


+ — 




D a 


M n 


n 


M t 


M 




M 



D n 



V0.2 



10- 3 



2/3 



(47) 



we see that condition (iii) is satisfied if 

1/2 



t D >7.3 x 10 



\GM) 
\0.2J \W- 3 ) 



4/3 



■2/3 



Mi 
M 



10/9 



(48) 



Since to is on the order of a thousand years or more for 
orbital frequencies comparable to WD g-modes, condition 
(iii) is always satisfied for WD/ WD or WD/NS binaries. On 
the other hand, the conditions (i)-(iv) are not all satisfied 
for NS/NS or NS/BH binaries. 

For very massive companions, the inequality of equation 



implies (for M t /M ~ M'/M) 
M' 



(48 1 may not hold. Using equation ( 10 \ for to, equation (48 1 

7 / 1( y Q 



M 



<5.5 x 10 



M 



4 tUJg V 

\0.2/ 

9/4 



V 10 4 km I 



10- 3 

9/4 



6/5 



(49) 



The above inequality implies our estimates are valid for any 



8 J. Fuller and D. Lai 



feasible companion except a super-massive black holej^Ve 
can also use this inequality to examine the inaccuracy of our 
estimate in the highest mass cases of Figure [6] Figure [6] was 
generated using the n = 3 mode parameters listed in Table 
1 for a WD of M = 0.6M Q and R = 8.97 x 10 3 km. Plugging 



in these parameters, equation (491 requires 
M' 



< 1.7 x 10 3 



(50) 



for our analytical estimates in Figure [5]to be accurate. This 
explains why the analytical estimates of Figure [6] are ac- 
curate when M' <J 1000A/ but diverge from the numerical 
results when M' £ 1000M. 

Given the maximum mode amplitude reached during a 
resonance, we can now estimate the temporal duration of 
the resonance. Letting a a — c a e _! " at , the mode amplitude 
evolution equation |7| becomes 

z . (51) 



2 iio a c a 



u 3 



Assuming that during the resonance, u) a — mf! — 0, the 
right-hand-side of equation (51 1 can be taken as a constant, 
we then have 



iM'WlmQo 

2ui a D 3 



(52) 



i.e., the mode amplitude grows linearly in time. Thus, the 
duration of the resonance is of order 



3.42 



M' 



-1/3 



Mt 
M 



5/9 



-13/9^,-2/3 



3\l/2 



= 3.50 x 10 13 



M'\- 1/3 /MA 5/g 



0.2. 



-13/9 



M 



10- 



M 



-2/3 /^ 3N l/2 



(53) 



Since the dynamical time (7? 3 /GM) 1//2 for typical WDs is on 
the order of one second, the resonance duration is typically 
an hour or longer. Note that the above estimate is formally 
valid only when [u — m£l(D ct +)] t ICS < 1, so that we can 
set Lu a — mQ w for the duration of the resonance. Using 
equations (|47[) and (|53[l, we find [co a — mQ(D a+ )]t ICS « 1 for 



typical parameters. Thus we should consider equation (531 



as an order-of-magnitude estimate only. Also, we can check 
that the GW energy loss during the resonance, AE gw ~ 



(MM'/2D ol )(t tes /tn 
our derivation of E, 



is much less than E a 



a ,max 



, justifying 
given by equation (381. Indeed, the 



above condition simplifies to equation |49| since in both cases 
it is the energy carried away by gravitational waves that is 
limiting the mode growth. 

We can use the same method to solve for the size of 
the fluctuations in mode amplitude after resonance. Due to 
the symmetry of the harmonic oscillator, the fluctuation in 
mode amplitude about the mean value after the resonance is 



1 Obviously, our estimate would not apply for a WD in a highly 
eccentric orbit around an intermediate mass black hole, which 
may form in dense clusters as described by Ivanov & Papaloizou 
(2007). 




Figure 7. The amplitude of a mode as a function of time near 
its resonance for different values of the damping coefficient 7 a . 
The curves have *y a = (black), 7 Q = 10 -9 (dark blue), y a = 
10~ 8 (yellow-green), -y a = 10~ 7 (light blue), -y a = 10~ 6 (red), 
7 Q = 10 -5 (green), and f a = 10~ 4 (purple). For this mode, we 
have set u) a = 0.1s -1 and Q a = 1 X 10~ 4 so that t rcs fa 10 5 s. 
Note that the damping term does not greatly affect the maximum 
mode amplitude except when 7 Q ^> 10 -4 , or when 7 a u; a f.res ^ 1. 
Also note that modes with larger values of -f a evolve on a much 
shorter time scale because their orbits decay quickly due to the 
conversion of orbital energy into heat via mode damping. 



identical to the zeroth-order estimate of the mode amplitude 



before resonance [see eq. (28 1], i.e., 



Aa i 



M'W tm Q a 



(54) 



These fluctuations occur with frequency mQ — ui a , since this 
is the difference in frequency between the eigenfrequency 
at which the WD is oscillating and the orbital forcing fre- 
quency. So, as the orbital frequency continues to increase af- 
ter the resonance, the amplitude of the fluctuations becomes 
smaller while the frequency of the amplitude oscillations be- 
comes higher. 



6 EFFECT OF MODE DAMPING 

The results in the previous two sections neglect mode energy 
damping in the WD. Since the duration of the resonance is 



much longer than the mode period [see equation (53 1], inter- 



nal mode damping could affect the energy transfer during 
resonance if the damping rate is sufficiently large. To ad- 
dress this issue, we incorporate a phenomenological damp- 
ing rate —j a uj a a a to the mode equation (JtJ to study how 
mode damping affects energy transfer during a resonance. 
Figure [7] shows the excitation of a mode through resonance 
for different values of y a - We see that, as expected, when 
the internal damping time is larger than the resonance du- 



ration (equation 53 1, the maximum mode energy achieved 
in a resonance is unaffected. 

G-modes in white dwarfs are damped primarily by ra- 
diative diffusion. For sufficiently large mode amplitudes, 
non-linear damping is also important (e.g., Kumar & Good- 
man 1996: see section for more discussion on this issue). 
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Wu (1998) presents estimates for the non-adiabatic radia- 
tive damping rates of WD g-modes in terms of uii = juj t . 
Extrapolating Wu's values to I = 2 modes for a white dwarf 



of temperature T — 8720 K, we find 7 



10" 



for modes 



near n = 1 and 7 ~ 10~ for high-order modes with n J> 20. 
So, while the maximum amplitude of low-order modes is 
completely unaffected by non-adiabatic effects, high-order 
modes will damp on time scales similar to the excitation 
time scale. Therefore these high-order modes will attain am- 
plitudes somewhat smaller than estimated in the previous 
section. This continual process of high-order mode damping 
may extract energy out of the orbit more efficiently than 
discrete resonance events, causing a steady decay of the bi- 
nary's orbit. 



7 DISCUSSION 

We have shown that during the orbital decay of compact 
white dwarf binaries (WD/WD, WD/NS or WD/BH), a se- 
ries of g-modes can be tidally excited to large amplitudes 
(up to 0.1 in dimensionlcss units) as the orbital frequency 
sweeps through the resonant mode frequencies. Such mode 
excitations can significantly affect the orbital decay rate near 
resonance. Indeed, to properly calculate the resonant mode 
amplitude, it is necessary to take into account of the back- 
reaction of the excited modes on the orbit. One consequence 
of the resonant mode excitations is that the low-frequency 
(^ 10 -2 Hz) gravitational waveforms emitted by the binary, 
detectable by LISA, will deviate significantly from the point- 
mass binary prediction. This is in contrast to the case of 
neutron star binaries (NS/NS or NS/BH) studied previously 
(Reisenegger & Goldreich 1994; Lai 1994; Shibata 1994; Ho 
& Lai 1999; Lai & Wu 2006; Flanagan & Racine 2006), where 
the resonant mode amplitude is normally too small to affect 
the binary orbital decay rate and the gravitational wave- 
forms to be detected by ground-based gravitational wave 
detectors such as LIGO and VIRGO. 

In the case of WD binaries studied in this paper, the 
number of of orbits skipped as a result of a resonant mode 
excitation is 

to E, 



AN orh = -^^^, (55) 

J orb -L'orb 

where to is the gravitational wave decay time scale given 



in equation (10 1, P or b is the orbital period, and i?orb is the 



orbital energy at resonance. Using equation (381 for 13 ajmax , 
we find 



A7V orb =3.4 x 10' 



(— ) 

\10 4 km/ 



( 

2/3 / x 2/3 

lfp 



1/9 



n 



OAs' 



■11/9 



(56) 



The number of skipped orbital cycles should be compared 
to the number of orbits in a decay time, expressed by 



dN 



1 



d In Q 3ir 



Qtr 



= 3.3 x 10 



s ( Afg \( M t 



\MM> J V 2M< 



1/3 



0.1s- 1 



-5/3 



(57) 



The huge number of skipped orbital cycles implies that such 
a resonant interaction would be important, but because the 



number of skipped orbital cycles is much smaller than the 
number of orbital cycles in a decay time, resonances will not 
dominate the decay process. 

A second possible consequence of resonant mode exci- 
tations is that the large mode energy may lead to significant 
heating of the white dwarf prior to the binary merger. In- 



deed, equation ( 38 1 shows that for typical binary parame- 
ters, the mode energy can be a significant fraction (~ 10 -4 - 
10 -3 ) of the gravitational binding energy of the star, and 
comparable to the thermal energy. Indeed, the thermal en- 
ergy of the WD is of order E t h ~ ^ , where T c is the core 
temperature of the WD and A is the mean atomic weight. 
The ratio of post- resonance mode energy to thermal energy 
is then 



E 



a . max 



Eth 



1.7 
Q 



~MJ \~M 

2/3 /1Q 7 K 



-8/9 



0.2/ 



10/9 



10- 



Tr. 



M_\ A0 4 km 



R 



(58) 



This implies that the white dwarf may become bright thou- 
sands of years before binary merger. 

A third consequence that may result from a resonance is 
significant spin-up of the WD. If we assume that all the an- 
gular momentum transferred to the WD during a resonance 
eventually manifests as rigid body rotation of the WD, the 
change in spin frequency of the WD is 



Ea- 



rn 



(59) 



where Q s is the spin frequency of the WD and / is its mo- 
ment of inertia. Plugging in our expression for E atUliix , we 
find 

,'0.2\ /M'\ 4/ VM^- 8/9 



(£f /9 ( 



V10- 



2/3 



n. 



(60) 



where n = I /(MR 2 ) w 0.2. We can thus see that a given 
resonance may deposit enough angular momentum to com- 
pletely spin up the WD (or significantly alter its spin) by 
the time the mode damps out. This implies that mode res- 
onances are potentially very important in the spin synchro- 
nization process. 

However, before these implications can be taken seri- 
ously, one should be aware of the limitations of the present 
study. One issue is the assumption that the white dwarf 
is non-rotating (and not synchronized), already commented 
on in section 1. More importantly, we have assumed that 
the white dwarf oscillations can be calculated in the linear 
regime. While the mass-averaged dimensionless amplitude 
|oa| = \b a \ of the excited g-mode is less than 0.1 [see Fig. 5 



and Eq. (39 1], the physical fluid displacement in the stel- 



lar envelope is much larger since g-modes of white dwarfs 
are mainly concentrated in the outer, non-degenerate layers. 
Figure [8] gives some examples: it shows the horizontal and 
radial displacements of three modes at their post-resonance 
amplitudes. These are obtained from £ = a a £ a , with £ a the 
normalized eigenfunction (see section |3| and a a computed 
from equation (391 with M' = M. In general, the linear ap- 
proximation is valid only if j£j <C |£v| -1 , where k r is the 
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0.0 0.2 0.4 0.6 0.8 1.0 
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Figure 8. The horizontal (solid line) and radial (dashed line) 
displacements of the n = 1 (top panel) and n = 5 (bottom panel) 
modes as a function of radius. The physical displacements are 
calculated using the analytical estimates for the post-resonance 
amplitudes given in equation |39| using M = M', giving |ai| = 
0.0566 and \a^\ = 0.0695. Also shown is inverse of the WKB wave 
number l/k r (dotted line). 

WKB wave number, given by 

fc 2 ^ + l)(iV 2 -^ 

Clearly, the three modes depicted in Fig. 8 strongly violate 
the linear approximation beyond the radius r ~ 0.85-R, near 
the jump in iV 2 associated with the carbon-helium bound- 
ary. 

Therefore, the results presented in this paper should 
be treated with caution as nonlinear effects will likely limit 
mode growth. Rather than increasing to the large displace- 
ments shown in Figure [8] the white dwarf oscillations will 
undergo non-linear processes such as mode coupling that 
will transfer energy to high-order modes. These high-order 
modes have much shorter wavelengths and thus damp on 
very short time scales. As nonlinearity is most important in 
the outer layers of the white dwarf, we expect that the ex- 
cited oscillation will dissipate its energy preferentially in this 
outer layer and will not reflect back into the stellar interior. 
We plan to address these issues in our next paper (Fuller & 
Lai 2010, in preparation). 
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